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1 Introduction 



The Calogero Model (CM) [1] - [6] is a well-known exactly solvable many-body sys- 
tem, both at the classical and quantum levels. The CM and its various descendants 
continue to draw considerable interest due to their diverse physical applications in 
systems such as random matrices [7], fractional statistics [8]- [11], gravity and black 
hole physics [12, 13], spin chains [14, 15], sohtons [16]-[20], 2D Yang-Mills theory 
[21, 22], lowest Landau level (LLL) anyon models [23, 24], Chern-Simons matrix 
model [25]-[27], Laughhn-Hall states [28]-[30], and unoriented superstrings in two 
dimensions [31]. 

Calogero 's original model describes N indistinguishable particles on the line 
which interact through an inverse-square two-body interaction. It is well-known, 
however, that the CM may alternatively be interpreted in terms of N free particles 
obeying generalized exclusion statistics [9, 10, 32, 33]. 

Haldane's formulation of statistics may be extended to systems made of different 
species of particles, in which the interspecies statistical coupling depends on the 
species being coupled. This may be implemented in a multi-species generalization 
of the CM in which particles have different masses and different couplings to each 
other [34]-[38]. 

Quite a few such generalized multi-sepcies Calogero models exist, but contrary 
to the original CM, knowledge about their exact solvabihty was rather tenuous. The 
recent breakthrough in this front derives from the papers [39]- [42]. The authors of 
[39] introduced deformed Calogero models, related to root systems of superalgebras, 
and gave effectively a proof of their integrability. In [40] they presented a more 



2 



conceptual proof by using shifted super- Jack polynomials. In related developments, 
the authors of [41, 42] introduced a supersymmetric generalizations of the CM which 
was based on Jacobians for the radial coordinates on certain superspaces. Both 
aforementioned models are closely related to the multi-family generalization of the 
CM introduced in [43, 44]. 

Motivated by these developments, in the present paper we investigate the latter 
model in the limit in which each family contains a large number of particles. In this 
hmit, the high-density hmit, the system is amenable to large- A?" collective-field for- 
mulation. As is well-known, the collective theory offers a continuum field-theoretic 
framework for studying interesting aspects of many-particle systems, somwhat anal- 
ogous to the continuum hydrodynamic description of fluids. It is appropriate to 
mention at this point the recent review on the collective-field and other continuum 
approaches to the spin-Calogero-Sutherland model [45]. The collective formulation 
has several virtues. In the large- A?" limit dynamics is governed by saddle points of 
the effective collective action, which contains the leading quantum effects. Thus, by 
extremizing this action, one is able to compute the uniform ground-state collective- 
field configuration, as well as topological and non-topological soliton configurations, 
and their corresponding energy eigenvalues. By expanding around these extrema 
it is possible to go beyond the large- A?^ leading order and obtain the spectrum and 
wave-functionals of the quadratice fiuctuations around these semiclassical configura- 
tions, and also to compute the corresponding density-density correlation functions. 

Beside these obvious advantages, the collective-field theory provides a natural 
framework for analyzing symmetries of the system which cannot be seen directly 
in the original (finite) A'"-particle quantum system. An important example in this 
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respect is the strong-weak coupling duality symmetry of the one-family Calogero 
model discussed in [46]. In this paper we generalize this approach to the multi- 
family Calogero model. We show that the resulting collective field theory is invari- 
ant under certain duality transformations, which interchange, among other things, 
particles and antiparticles, and thus generalize the duality symmetry [46] of the 
ordinary Calogero model. We identify all these duahties, which form an Abelian 
group, and study their consequences. In particular, the investigations carried in this 
paper will enable us to find the conditions under which collective quasi-particles 
describing density fluctuations in the the F-family Calogero model can be identified 
with those of an effective one-family Calogero model. As a by-product, this may 
help to better understand the exact solvability of some of the recently proposed 
two- family Calogero models [3 9] -[42]. We stress that the duality relations derived 
and discussed in this paper are exact symmetries of the collective-field hamiltonian, 
as opposed to the approxiamte duality symmetries discussed in [60, 61] . 

This paper is organized as follows: In Section 2 we review some known facts 
about the two-family Calogero model [47]. We then derive the quantum collective- 
field Hamiltonian, including divergent terms which are inherent to the collective-field 
formulation. In Section 3 we show that this collective field theory is invariant under 
certain duality transformations, which interchange, among other things, particles 
and antiparticles, and thus generalize the well-known strong- weak coupling duality 
symmetry of the ordinary Calogero model. We make a complete list of these du- 
alities, show that they form an Abelian group, and study their consequences. In 
Section 4 we concentrate on the leading large- A?" behavior of the collective Hamilto- 
nian of the two-family model and identify the effective potential which dominates 
dynamics in this limit. The large- A^ uniform ground state configuration is then 
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found by minimizing the effective potential, and the fluctuations around it are then 
investigated. The Hamiltonian which governs these small quadratic density fluctua- 
tions contains divergent terms. By carefully choosing the creation and annihilation 
operators describing these excitations, these divergent terms can be renormalized 
away. Section 5 is devoted to diagonalizing the latter Hamiltonian, resulting in two 
c?ecoup/ed quasi-particle dispersions. One of these dispersions, at small momenta, has 
the structure of a single-species Calogero model with some effective particle density, 
mass and statistical parameter, which we calculate. The other dispersion is that of 
free massive particles. In Section 6 we consider the ground-state wave-functional of 
small fluctuations, and calculate the corresponding static density-density correlation 
functions of the two-family model. Not surprisingly, the long-distance behavior of 
these correlation functions coincides with that of the effective single-species Calogero 
model alluded to above. In Section 7, we generalize our study of the two-family case 
to the F-family Calogero model (with no three-body interactions). We find all the 
exact duality transformations which leave the collective Hamiltonian of the F-family 
Calogero model invariant, and identify the Abelian group they form. We also obtain 
the ground-state and study the small fluctuations around it. In particular, we com- 
pute the low-momentum behavior of the dispersion relations of the quasi-particles 
corresponding to these fluctuations. As in the two-family case, one of these dis- 
persions has the structure of a single-species Calogero model with some effective 
particle density, mass and statistical parameter. The remaining F — 1 dispersions 
are those of free massive particles. We then compare our dispersion laws with those 
found before in the generalized Thomas- Fermi approach [36]. Finally, in Section 8 
we argue how our collective-field theory approach may explain the exact solvability 
of the two-family Calogero models proposed in recent papers. 
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2 Two-family collective-field Hamiltonian 



In order to get oriented, and in preparation for studying the multi-species Calogero 
model, let us start by recalling the essential features of the two-family Calogero 
system. The Hamiltonian of this model reads [44] 

Here, the first family contains A^i particles of mass m\ at positions Xi^ i — 
l,2,...,Ni, and the second one contains N2 particles of mass 777.2 at positions 
Xa, OL = 1,2,...,A^2- AH particles interact via two-body inverse-square potentials. 
The interaction strengths within each family are parametrized by the coupling con- 
stants Ai and A2, respectively. The interaction strength between particles of the 
first and the second family is parametrized by A12. For reasons that should become 
clear in later sections, we shall assume in this paper that Ai and A2 are always posi- 
tive, whereas the sign of A12 is unrestricted. In any case, these couplings are always 
assumed to be set such that (1) has a well-defined ground state. 

Note that we did not include in (1) a confining potential. This is not really a 
problem, as we can always add a very shallow confining potential to regulate the 
problem (in the case of purely repulsive interactions), or else, consider the particles 
confined to a very large circle (i.e., consider (1) as the large radius limit of the 
Calogero-Sutherland model [2]). We shall henceforth tacitly assume that the system 
is thus properly regularized at large distances. 

In (1) we imposed the restriction that there be no three-body interactions, which 
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requires [43, 44, 48, 49, 50] 

Ai A2 A12 



mi^ 17111712 



(2) 



This particular relation was also displayed in [36], as a consequence of the require- 
ment that the asymptotic Bethe Ansatz should be applicable to the ground state of 
the multi-species model. It follows from (2) that 

^^12 — '^l'^2 ■ (3) 

We assume that (3) holds throughout this paper. The Hamiltonian (1) describes 
the simplest multi-species Calogero model for particles on the line, interacting only 
with two-body potentials. The singularities of the Hamiltonian (1) at points where 
particles coincide implies that the many-body eigenfunctions contain a product of 
Jastrow-type pref actors 1111121112, where 

Ni 

Hi = n - Xj)^^ , 

i<j 

N2 

n2 = l[i^--xfi)'\ 

Ni,N2 

ni2 = n (^.-^a)'"- (4) 
i,a 

These Jastrow factors vanish (for positive A's) at particle coincidence points, and 
multiply that part of the wave-function which is totally symmetric under any per- 
mutation of identical particles. It is precisely these symmetric wave-functions on 
which the collective field operators act, as explained below. 

Let us make a couple of elementary comments at this point: First, we note from 
(4) that for A^ = and A^ = 1, the model describes two families of interacting 
bosons and fermions, respectively. Second, superficial glance at (1) may suggest 



that the model breaks into two independent Calogero systems at A12 = 1 . This 
is wrong, however. The two famihes remain correlated at A12 = 1, and display 
mutual fermionic behavior. Namely, the relevant prefactor in (4) picks up a factor 
(—1)^^^ under the exchange of particle indices {i ^ a). We will see, after making 
the trasformation to collective variables, that this gives a non-trivial interacting 
bosonic theory of the two collective fields describing different families of particles. 

Let us recall at this point some of the basic ideas of the collective-field method 
(adapted for our two- family Calogero model [47]): Instead of solving the Schrodinger 
equation associated with (1) for the many-body eigenf unctions, subjected to the ap- 
propriate particle statistics (Bosonic, Fermionic of fractional), we restrict ourselves 
to functions which are totally symmetric under any permutation of identical parti- 
cles. This we achieve by stripping off the Jastrow factors (4) from the eigenf unctions, 
which means performing on (1) the similarity transformation 

H^H = nr2^n2-^nr^i/nin2ni2 , (5) 

where the Hamiltonian 

2mi ^ dxi"^ mi y ^ Xi - Xj j;^ - Xa j dxi 

I N2 Q2 I f I 1 \ d 

-^E/t- — ^2E^^ + Ai2E— (6) 

2rn2 ^1 OXa^ m2 y e^/s^a-Xp - Xi J dx^ 

Note that H does not contain the singular two-body interactions. By construction, 
this Hamiltonian is hermitian with respect to the measure 

diJ,{xi, Xa) = (ninallia)^ d^^Xi d^^Xa , 

(as opposed to the original Hamiltonian H in (1), which is hermitian with respect 
to the fiat Cartesian measure) . 



We can think of the symmetric many-body wave-functions acted upon by H as 
functions depending on all possible symmetric combinations of particle coordinates 
(symmetric separately in the coordinates of each family of particles) . These combi- 
nations form an ovcrcomplete set of variables. However, as explained below, in the 
continuum limit, redundancy of these symmetric variables has a negligible effect. 
The set of these symmetric variables can be generated, for example, by producs of 
moments of the collective - or density - fields 

Ni N2 

Pl(^) = ~ ^i) ' P2{x) ^ ^ S{X - Xa) . (7) 

i=l a=l 

The collective-field theory for the two-family Calogero model is obtained by chang- 
ing variables from the particle coordinates Xi and Xa to the density fields pi (x) 
and P2{x) . This transformation replaces the finitely many variables Xj and x^ 
by two continuous fields, which is just another manifestation of overcompleteness 
of the collective variables. Clearly, description of the particle systems in terms of 
continuous fields becomes an effectively good description in the high density limit. 
In this limit the mean interparticle distance is much smaller than than any relevant 
physical length-scale, and the S- spikes in (7) can be smoothed out into well-behaved 
countinuum fields. All this is in direct analogy to the hydrodynamical effective de- 
scription of fluids, which replaces the microscopic atomistic formulation. Of course, 
the large density limit means that we have taken the large- A^i,A^2 hmit. (It is 
understood throughout this paper that A^i and N2 tend to infinity at comparable 
rates.) 

Changing variables from particle coordinates Xi,Xa to the collective fields (7) 
implies that we should express all partial derivatives in the Hamiltonian H in (6) as 

— ^[dx^^-^, -^=[dx^^-^, (8) 
dxi J dxi Spi{x) ' dxa J dxa Sp2{x) ' 
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where we applied the differentiation chain rule. 



In the large — A'^i, A'2 limit, the Hamiltonian H can be expressed entirely in terms 
of the collective fields pi, p2 and their canonical conjugate momenta 

T^li.^) = ' ^2(3;) = -^7-^ , (9) 

dpi[x) dp2[x) 
as we show below. It follows from (8) and (9) that the particle momentum operators 

(acting on symmetric wave-functions) may be expressed in terms of the collective- 
field momenta at particular points on the line as 

Pi = -^[{Xi) , Pa = -A{Xa) , (10) 

(where n'^{x) = d^Traix)). Finally, note from (7) that the collective fields obey the 
normalization conditions 

J dxpi{x) = Ni, J dxp2{x) = N2 . (11) 

The density fields Pi, P2, and their conjugate momenta tti,tt2, satisfy the 
equal-time canonical commutation relations^ 

[pa{x),7rb{y)] = iSabS{x -y) , a, 6 =1,2, (12) 

[Pa{x),pb{y)] = ['Ka{x),7rb{y)] =0. 

By substituting (7)- (10) in (6), we obtain the continuum-limit expression for H 

as 



-if 

H ^ 2rn~ J ^^'^i(^)(^^^i(^))' 



^According to (11), the zero-momentum modes of the density fields are constrained, i.e., non- 
dynamical. This affects the first set of commutation relations in (12), whose precise form is 
[pa{x) , Trb{y)] = iSab{S{x — y) — (l/0)> where I is the size of the large one-dimensional box in which 
the system is quantized, which is much larger than the macroscopic size L of the particle condensate 
in the system. In what follows, we can safely ignore this l/l correction in the commutation relations. 
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i f, I \( ^^ -^ d^Pi , , f dypijy) , , f dyp2{y) \ ^ 

/ dxpi{x) — h Aif \-Xi2f d^ni{x) 

mi J \2 J ^ — y J X — y J 

+ 2^/ dxp2{x){d^T:2{x)f 
^ r \ (^^-'^9xP2 . . J'dyp2{y) fdypi{y)\^ .,„^ 

/ dxp2{x) h A2f h A12-/- 9a;7r2(x), (13) 

m2J \2 P2 J X - y J x-y I 

where ^ denotes Cauchy's principal value 



It can be shown [51] that (13) is hermitian with respect to the functional measure^ 
T^p[pu P2] = J\pi, P2] n dpi{x) dp2{x) , (14) 

X 

where J[pi, P2] is the Jacobian of the transformation from {xi, Xa} to the collective 
fields {pi{x) , p2{x)} . In the large - Ni,N2 limit it is given by [47] 

lnJ=(l-Ai) j dxpi{x)\npi{x) + {1 - X2) j dxp2{x)\n p2{x) 

-2Ai2 j dxdy pi{x) In \x - y\p2{y). (15) 

It is more convenient to work with an Hamiltonian, which unlike (13), is hermitian 
with respect to the fiat functional Cartesian measure Yix dpi{x) dp2{x) . This we 
achieve by means of the similarity transformation ijj — > J^i/j , H — > Hcou = J^HJ~^ , 
where the continuum collective Hamiltonian is 

Hcoll 

'Xi-ld^Pi fdypiiy) , , fdyp2{yy 



2m 

+ 



1 f, //^ / ^//^, 1 /■> / jXi-ld^pi ..fdypiiy) fdyp2{ 

— / dx-K-^^x) pi{x)t:^{x)+- — / dxpi[x)\ — h Aif h Ai2f 

Til J 2mi J \2 pi J X — y J x — 



- J dxn,{x) P2{X) n2{x) + ^ J dxp2{x) [-^—- + + Al2/^ 



2m2 J 2m2 J \2 p2 J x — y J x — y 



^By definition (recall (7)), this measure is defined only over positive values of pi, P2 • 



where Hsmg denotes a singular boundary contribution^ : 

Hsina — — ( dx \ -^^Plix) + ^^/Oofx) I 8^ 

''"^ J \2mx'^^ ' 2m2 ^ 7 x-y 



(16) 



y=x 



and P is the principal part symbol. 

It is well-known [52] that to leading order in the expansion, collective dynam- 
ics of our system (properly represented by well-behaved observables) is determined 
by the classical equations of motion resulting from (16). 

As a concluding remark, note that the terms in (16) are independent of 

the normalization (11) of pi and p2. Thus, they are clearly of next-to-leading order 
in the expansion relative to the other terms. Nevertheless, these terms play 
an important role in in dynamics of the model, e.g., in governing small fluctuations 
around the ground state (see Sections 4 and 5), in the soliton sector of the model 
[47], and in establishing the duality symmetries of (16), which we do in the next 
section. 

3 Duality transformations and symmetries 

In this section we introduce and study the duality transformation and the cor- 
responding symmetries of the collective- field Hamiltonian (16). It is straightfor- 



dx — 

X — y 



Note that the singular coefficients of pi and p2 in (17) are in fact independent of x, as evidently, 
dx ^\x=o ~ ' ^^6^6 the last step we used the standard representation ^ = 



J/=X 



x^+e^ 
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ward to check that the Hamiltonian (16) is invariant under the following set of 
transformations^ of the parameters: 

r 1 r 1 ~ mi ^ ma r 1 , . 

Ai = — , A2 = — , mi = — — , 1712 = — 7-, A12 = — (18) 
Ai A2 Ai A2 A12 

and of the operators: 

Pi = -AlPl, P2 = -A2P2, TTl = -T^, 7^2 = -T^ ■ (19) 

Al A2 

Derivation of these transformations, as well as of the other sets of transformations 
discussed below in this section is straightforward: Consider two two-family Calogero 
systems with collective Hamihonians of the form (16). Assume that the tildcd quan- 
tities, corresponding say, to the second system, are simple homogeneous functions of 
the untilded ones (allowing sign flips). Assume further that the tilded and untilded 
quantities appearing in (2) have the same common value. The transformations (18) 
and (19) then follow as the unique solution of equating the two Hamiltonians, which 
is neither an identity, nor a trivial permutation. (The (Aj — l)dxPi/ Pi terms in (16), 
are crucial in obtaining these transformations uniquely.) 

Let us denote the set of transformations (18) and (19) by T12. These transforma- 
tions are canonical, as they preserve the commutation relations (12). For obvious 
reasons, we refer to the transformations T12 as the strong-weak coupling duality 
transformation. Thus, we see that our Hamiltonian, expressed in terms of the new 
tilded parameters and operators, is identical in form to the original one, but with 
Al and A2 and the inter-family coupling A12 turned into their reciprocal values; 
with A^i and A^2 turned, respectively, into N\ — —XiNi and N2 — —X2N2 
and, finally, with masses mi and 1712 turned into — ^ and — The minus 



*Note that (18) and (19) do not constitute a symmetry of the original Hamiltonian (1). 
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signs which occur in these identifications are all important: By drawing analogy to 
a similar situation in the one-family case [46, 53], we interpret all negative values 
of the parameters and densities as those pertaining to holes, or anti-particles. Now, 
strictly speaking, since iVj and are integers, this interpretation is consistent only 
for rational couplings, as was discussed in [46, 53]. 

We further note that (16) is invariant also under two more sets of canonical 
duality transformations. The first one, which we denote by Ti, is comprised of the 
set of transformations of parameters 

Ai = Ai, A2 = ^, 77li = mi, 7712 = -^, A12 = (20) 

A2 A2 A2 

and of the operators 

Pi = Pi, P2 = -A2P2, TTl = TTl, 7f2 = -— . (21) 

A2 

Negative values of masses, densities and momenta, as in the previous case, refer 
to holes. These transformations map the two-family Calogero model of particles 
(positive 1711,1712, Pi and p2 ) with inter-family interaction strength A12 into the 
dual two-family Calogero model of particles {mi, pi) and holes {rh2, P2) with 
the inter-family interaction strength — The second (and last) set of duahty 
symmetries of (16), which we denote by T2, is obtained from (20) and (21) simply 
by permuting the family indices 1 -e^ 2 . 

It is easy to check that the duality transformations Ti, T2, T12, together with 
the identity transformation I, form an Abelian group under composition, in which 
each element squares to I, and where T1T2 = T12, T1T12 = T2 and T2T12 = Ti. This 
is readily identified as Klein's four-group. The latter is isomorphic to Z2 Z2 , 
where the two Z2 factors are {/, Ti} and {/, T2} . As we shall see in Section 7, this 
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group-theoretic pattern of duahties persists also in the generic case. 

3.1 Resemblence of special two- family Calogero models to 
single-family Calogero models 

As an interesting application of the duality symmetry group, consider the special 
case of the two-family Calogero model (1) in which^ ^2 = and m2 = ■ From 
(2) we then find that A12 — — —1. Note that the two particle families in 

this system are generically manifestly distinct. Nevertheless, this distinction is, in 
some sense, an illusion. To see this, note that by the duality transformations (20) 
and (21), namely, the element Ti of the duality group, this system is equivalent to 
a two-family system with parameters Ai = A2 = A12 = Ai , rhi — 7712 = mi and 
densities pi = pi , p2 = ■ In the latter dual system, the two famihes are identical! 
For this reason, we may refer to the two dimensional locus 

A2 = ^, m2 = -^, Ai2 = -1 (22) 

in parameter space as the "surface of hidden identity" , or SOHI. Thus, the special 
two-family Calogero model we started with resembles the single-family Calogero 
model specified by 

A = Al , m = mi , p = Pi + P2 = Pi - ^P2 • (23) 

Al 

Similarly, by inverting the roles of family indices 1 -s-^ 2 in the previous discussion, 
which leaves us on the SOHI (22), and then applying the duality tranformation 

T2, we shall conclude that the special two-family Calogero model we started with 

^If we ignore inter-family coupling, we can think of this system as made of two single-family 
models, related by the one-family version of the strong-weak coupling duality, save for the relation 
between densities. 
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resembles the single-family Calogero model specified by 

A = A2 , 171 = 7712, P = P2- ^Pl ■ (24) 

In both cases, the effective single-species collective field p actually shares the 
statistics A and the mass m with the first or the second family, respectively. Note 
that these two cases can be mapped one each other by the duality transformation 
T12. Thus, the SOHl (22) is left invariant under T12. However, the latter does not 
act on it freely, as Ai = A2 = — A12 = 1 and mi = — m2 is a fixed line. Models lying 
on this line are comprised of particles and their antiparticles, and only particles and 
antiparticlcs interact (repulsively). 

Note that we described the relation between the original special two-family mod- 
els and the corresponding single-family models merely as "resemblance" . They are 
certainly not identical! The density operator p appearing on the LHS of (23) and 
(24), which corresponds to the single- family Calogero model, is defined in a Hilbert 
space made of many-body wave functions which are completely symmetric in the 
coordinates of all particles, pi and p2, on the other hand, are symmetric only in the 
coordinates of particles of each family separately. The best one could do is perhaps 
to consider the two-family system with identical families (the one dual to the special 
two-family systems we started with) as a one-family system divided into two parts, 
differing by some internal quantum number, in which one symmetrizes in each sector 
separately. However, this means one should also contrapt an actual physical context 
to justify such separate symmetrization. 
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4 Ground state and small fluctuations 



In this section we shall concentrate on the ground-state of the two-family Calogero 
system and the small excitations above it. To this end, it is sufficient to consider 
only static density fields, for which, of course, tti = 7T2 = 0. Evaluating the Hamilto- 
nian (16) on such fields and imposing the normalization contraints (11), yields the 
effective potential, given (to leading order in the 1/Ni and I/N2 expansion) by 
1 /■ , / ./Ai-l(9^pi fdypi{y) fdyp2{y)\ 



V 



zrrii J \ ^ Pi J X — y J x — y J 



2m2 J \2 p2 J x-y J x-y J 

+pi (^Ni - J dxpi{x)^ + p2 (^N2 - J dxp2{x)j . (25) 



Here pi and p2 are the chemical potentials which impose (11).^ 

The expression in (25) contains bilocal as well as trilocal terms in the densities. 

We can avoid the menacing trilocal terms by applying standard tricks as follows: As 

was mentioned at the beginning of the previous section, we tacitly assume that the 

system is properly regularized at large distances, meaning in particular, that pi{x) 

and p2{x) are always of compact support. The principal value distribution, acting 

on such functions, satisfies the identity 

P P P P P P 

— + — — + — — = n^S(x - y)S(x - z) . (26) 

x—yx—z y—zy—x z—xz—y 

Making use of (26) in (25), we obtain 



V = — / dx{Xipi + A12P2) + 

oAimi J 



^In (25) we have absorbed the constant singular coefficients of pi and p2 which appear in Hging 
in (17) as infinite shifts in fii and /X2, and also omitted the infinite, density independent constant 

in Hsing- 
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zmi J J X — y 2m2 J J x — y 

J J X — y zrrio J J x — 



2mi i J X — y 2m2 J J x — y 

+Pi (iVi - I dxpi{x)^ + 112 {n2 - I dxp2{x)^ . (27) 

This expression for V is evidently devoid of trilocal terms. In Appendix A we present 
an alternative derivation of (27). 

As was mentioned at the end of the previous section, in the large- A'" limit, dy- 
namics of the system is semiclassical, governed by the saddle points of (16). In 
particular, the ground-state energy and the corresponding density distributions can 
be found, to leading order in by varying the effective potential (27) with respect 
to pi and p2, which yields the two coupled integro-differential, non-linear equations 

.2 (Ai - if d^pi {Xi - if d^cPi^ 
{MP1 + M2P2) ( ) ) 



2mi 8mi pi Ami Pi 

Ai(Ai - 1) 



[Idy^-lEM _ bl^h^ + h^)Uy^-yPM ^ (28) 
J x — V 2 ■m^ 1712 J X — y 



mi J X — y 2 mi m2 J x — y 

7r%, , , x2 {X2-lf.d^p2f {X2-lf r. ^d^p2. 
-(A2P2 + A12P1) — (— ) ' 



2m2 8m2 P2 4m2 P2 

_ X2{\2 - I) f^ydypM _ A12 (Ai_^ ^ ^)irf,^PlM ^ (29) 
m2 X — y 2 mi m2 J x — y 

In the ground state, the two particle types evidently condense into uniform density 



configurations 



^1 ^2 
Pi,o = , P2,o = , (30) 



inside the spatial box of length L, which confines the particles. These uniform 
densities are clearly solutions of the variational equations (28) and (29), since for 
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uniform densities the latter equations reduce simply to relations which determine 
the chemical potentials. To conclude our discussion of the semiclassical limit, note 
that to leading order in the ground-state energy corresponding to (30) follows 
from (27): 

This is, in fact, the exact result for finite A^i and N2 [36]. It is worth mentioning 

here that (28) and (29) possess also nonuniform static soliton solutions, which were 
studied extensively in [47, 54]. 

The collective-field formalism provides a systematic framework for the ex- 
pansion [55, 56]. In particular, by expanding around the classical uniform configu- 
rations (30) (with their corresponding null momenta), we can go beyond the leading 
order and obtain the spectrum of low-lying excitations above the ground state. Here 
we shall concentrate on the next-to-leading terms. To this end we write 

pi{x) = pi^o + , P2{x) ^ p2,o + mix) , (32) 

where rii{x) and i]2{x) arc small density fluctuations, with rja being typically of order 
l/Na- Substituting (32) in the Hamiltonian (16), and expanding it to quadratic order 
in TTi, 7r2, r^i , and 772 , we obtain 

2m2 J 2m2 J \ 2 p2,o J x — y J x — y J 

+Hsing + Eq . (33) 

The quadratic approximation (33) to H^ou, as the original collective Hamiltonian 
(16), is invariant under the duality transformations Ti, T^, and T12 in Section 3. Eq. 
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(33) governs the leading-order fluctuations of the system around the ground state. 
Its quadratic structure clearly calls for introducing the two operators 

Ai[x) = d^'Ki[x)+i — h Aif h Ai2f , (34) 

V 2 pi,o J x-y J x-y J 

Af\ . [ ^2- l d^r]2 , ^ f dymjy) , , f dyr]i{y) \ 
A2{x) = dx'K2{x) + i — h Asf h \i2f , (35) 

V 2 p2,o J x-y J x-y J 

along with their hermitian adjoints. Evidently, these operators are the key to com- 
puting the collective quasiparticle spectrum of small fluctuations around the ground 
state. Higher order terms in the expansion of Hcoii around the uniform densities 
yield the interactions among these quasiparticles. 

It can be checked that the duality transformations of Section 3 act on Ai{x) and 
A2{x) as follows: Ai is invariant under Ti while T1A2 — > —j^A2 ; A2 is invariant un- 
der T2 while T2A1 — > —^7^1 ; flnally, ^12^1,2 — ~ a7I^i,2 (the latter transformation 
is consistent with the group multiplication T1T2 = T12 ) . 

The operators (34), (35) and their adjoints satisfy the commutation relations 



1 - Ala o ^/ N ^^ a P 



[A,{x), Al{y)] = -d,dy6{x -y) + 2\^d, 



Pi,o ^ ' ■ ""x-y' 

[A2{x), Al{y)] = ^i^d,dyS{x -y) + 2X28,^, (36) 
P2,o X-y 

[A2{x), A\{y)] = [A,{x), 4(y)] = 2Ai2a ^ 



X -I 

x-y 



with all other commutators vanishing. In terms of these operators the quadratic 
Hamiltonian (33) takes the form 

Hcou = ^0 + £^ / dxA\{x)A,{x) + §^J dxAl{x)A2{x) 

dx[A,{x),A\{x)] + ^ / dx[A2{x),Al{x)]+K 
4m2 J 



4m I 



sing • 
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As can be seen from (36), the commutator terms in the last equation, which are 
essentially the sum over zero-point energies originating from the quadratic pieces in 
(33), cancel the the singular term Hging defined in (17). This is a typical behavior 
in collective field theory. Thus, we obtain the finite result 

Hcoii ^Eo + ^ f dxA{{x)A^{x) + ^ / dxAl{x)A2{x). (37) 
Zuii J 2fm2 J 

As our next step towards diagonalization of (37), we move to momentum space and 

express the operators A\ and A^ in terms of of their Fourier transforms 



A^{x) = f dke'^^a^ik) , A^ix) = f dke^'^a^ik) . (38) 

By definition of the duality trasnformations in Section 3, and from their action on 
Ai and A2 as described following (35), we obtain that their action on ai and 02 
reduces to either flipping a sign or not: oi is invariant under Ti while T\a2 —>■ —02 ; 
a2 is invariant under T2 while 720i — > — Oi ; finally, Ti2ai,2 ■ 



^From (34) and (35), and from the identity 

ikx 

-00 X — y 

it is easy to verify that 



/cxD pikx 
dx = ivre^'^^'signA; (39) 
-00 X — V 



ai{k) 



a2{k) 



/7rpi,o 
mi 



1712 



ikni{k) + ^- — —kfji{k) + n {Xirii{k) + Xi2fi2{k)) sign/c 



1 - A2 

ikTT2{k) + — kfj2{k) + n {X2'fi2{k) + Xi2fii{k)) signk 



(40) 



The Fourier modes 



Uk) = / ^e-'''Va(x) , Uk) = / ^^-'"^^aix) (41) 
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obey the commutation relations^'^ 

[Uk),MQ)]-^S{k + q), (42) 

with all other commutators being null. For convenience, we also write down the 
hermitian adjoints of (40): 



a\{k) = 



al{k) 



nil 



/ 7rp2,o 

1712 



-ikTTi{-k) + ^- — —kfii{-k) + TT {Xifji{-k) + Ai2?72(-^)) sign/c 



1 - A2 

-ikTT2{-k) + — kfj2{-k) + TT {X2fi2{-k) + Ai2?7i(-/c)) sign/c 



(43) 

Note that by combining the equations for ai,{k) and a\{—k) we can invert (40) and 
(43), and express r]i){k) and i^bik) as linear combinations of ai){k) ± a\{—k) . 

Finally, it follows from (40), (42) and (43) that*^ 

[aAk), a{ik')\ = { ^7^A;2 + ^i^|A';|^ Sik - k') = Ui(k)5(k - k'), (44) 
\ 2mi m-i J 

[a2(k), al(k')] = ( ^^^e + ^^^|A;|) 6{k - k') = u;2(k)5(k - k'), (45) 

V 2m2 7712 J 



[a,{k),al{k')] = [a2{k),a\{k')] = Xu^J ^^^\k\S{k - k') = uji2{k)5{k - k'), (46) 

y mim2 

and all other commutators vanish. For later use, it will be convenient to lump (44) 
- (46) into a real symmetric matrix [aa(/c), CL\{k')] — Uab{k)8{k — k') , with 

(47) 

y Ul2[k) UJ2[k) I 



^Duc to the fact that J dxrjaix) = 0, the zero-momentum modes of the fluctuation fields 

are nondynamical commuting numbers. See the footnote just above (12) for a similar matter 
concerning the full density fields. As was commented upon there, in the limit of very large spatial 

box, we can safely use (42). 

^Note that fil{k) = fja{—k) and T^aik) = Tra{—k), since r]a{x) and TTa{x) are hermitian field 
operators. 

^Alternatively, Eqs. (44) - (46) can be verified by Fourier-transforming (36) directly. 
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It is easy to verify that the duahty transformations of Section 3 act on uj as follows: 
Ti2 leaves it invariant, while both T\ and T2 transform it to a-iiou-}, (i.e., they flip 
the sign of the off-diagonal terms oj\2). All this is of course consistent with the 
action of these duality transformations on a\ and 02, as described following (38). 
Note in passing that we can infer a special case of the exact duality equivalence of 
systems with parameters on the SOHI (22) and two-family models made of identical 
families, as discussed in Section 3.1. This we achieve simply by demanding that 
<^ii(^) — ^22{k), which immediately leads to the pair of equations — and 



to begin with), wc obtain that the conditions defining the SOHI (22) comprise a 
solution of the first equation, which upon substitution into the second equation 
leads to the further relation p2,o = —^iPi,o (which did not arise in Section 3.1). 

In terms of the operators (40) and (43) the Hamiltonian (37) becomes 



In Section 5 we will show that the ground state of the fluctuation fields is annihiliated 
by ai{k) and 02 (/c). Thus, the ground-state energy Eq is not affected by the next- 
to-leading-order corrections stemming from the quantum collective-field fluctuations 
771 and 772. 

5 Diagonalization and dispersion laws 

Expression (48) is not really the diagonal form of the Hamiltonian, since the oper- 
ators ai and 02 do not commute - see (46). As usual, complete decomposition of 
(48) into independent normal modes can be achieved by properly rotating ai{k) and 



Aipi,o 
mi 



^ . Discarding the trivial solution (i.e., that the two families are identical 




(48) 
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a2{k) into two new commuting bosonic operators bi{k) and b2{k). Thus, we write 

ba{k) = j2uab{k)ab{k), (49) 

b=l 

where Uab{k) are c-number coefficients. They are chosen such that the matrix u 
diagonaiizes the commutator matrix (47). Since lu is reai and symmetric, u is reai 
orthogonal. Thus, 

viun^ ^n, (50) 



where 



Recall that the duahty transformations Ti, T2 and T12 transform u; either to itself 
or to (T^ucr^ , which is similar to co . Thus, the the eigenvalues ^aik) of (47) are 
invariant under these duality transformations. These eigenvalues are given by the 
standard formula 

ni,2{k) = ^ (u,{k) + U2{k) ± ^{ujr{k)-uj2{k)f + ^uj^2^{k)^ . (52) 

By construction, therefore, the new operators satisfy the decoupled commutation 
algebra 

[hi{k),h\{k')] = Qi{k)6{k-k'), 

[b2{k),bUk')]^Q2{k)S{k-k'), (53) 
h{k),bl{k')]^0. 

The coefficients Uab{k) of the orthogonal transformation can be expressed in terms 
of the rotation angle (j){k) as 

Uu = 1*22 = COS (p[k) = 
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Ui2 = -1*21 = sm (p[k) = = - (54j 

Note from (44) - (46) that the Uab{k) are all even functions of k, in addition to being 
real. A useful and simple identity which we mention in passing is 

cot 20 =5^^^. (55) 

By rescaling 

h{k) ^ ^J^hik) bi{k) , b2{k) ^ ^fihik) 62(A)), (56) 

we finally obtain the diagonal Hamiltonian 

H,M = ^0 + / dkn^{k)b\{k)bi{k) + I dkn2{k)bl{k)b2{k) . (57) 

fli{k) and fl2{k) are therefore the energy spectra of the physical fluctuations, i.e. 
the two decoupled collective oscillator modes, or quasiparticles of the system. 

As can be readily seen from the commutation algebra (44)- (46), the dispersions 
u!i{k) and u!2{k), evaluated at A12 — , actually represent the low-energy dispersions 
for the two non-interacting one-family Calogero systems. Stability of these latter 
two systems requires that the coefficients of the terms in uJi{k) and uj2{k) which 
are proportional to the leading term in the low- momentum behavior of these 
systems, be positive. Thus, we shall henceforth make the assumption that these 
terms are positive. Note further that these terms are invariant under the duality 
transformations of Section 3, and so are their signs. 

Since we know from [57] that the low-energy one-family dispersions u!i{k) and 
u!2{k) are correctly given only up to k^ (near k — 0), the same must be also true for 
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Qi{k) and (/c). By expanding the square root in (52) up to order k"^, we obtain^° 
(l-Ai)Aipi,o/mi^ + (l-A2)A2P2,oM2^/s^ , /Aipi,o , A2P2,o\|,i 

iZilfc) — TT+'^l 1 \ \k\ : 

Mpi,o/mi + X2P2,o/m2 2 \ rrii m2 J 

^jj^. ^ (1 - A2)Aipi,o/(mim2) + (1 - Ai)A2P2,o/(rnim2) k'^ .^gs 
^ Aipi,o/mi + A2P2,oM2 2 

The first dispersion fli{k) has the same structure as that of a single-species 
Calogero model [56] with some effective particle density p, mass m and in- 
teraction parameter A satisfying the following two relations: 

(1 - Ai)Aipi,o/mi2 + (1 - A2)A2P2,o/'^2^ 1 - A 



AiPi,o/ mi + A2P2,o/ m2 m 



(59) 



^lPl,0 ^ ^2P2,0 ^ Ap 
mi 7712 

The second dispersion Q,2{k) describes a free-particle. The effective mass m* of 
this quasi-particle is given by 

(1 - A2)Aipi,o/(mim2) + (1 - Ai)A2P2,o/(mim2) _ J_ 
Aipi,o/mi A2P2,o/^2 ' 

As was mentioned earlier, the exact expressions for Qi,2(A;) are invariant under the 
duality transformations of Section 3. It can be further checked that the approximate 
expressions in (58) retain this invariance. Thus, in particular, the expressions on the 
LHS's of (59) - (61) are all invariant under the duality transformations of Section 3. 

It is interesting to observe that for the special values of uniform densities con- 
figurations p2,o = Aipi^o; parameters A1A2 — A12 = 1, and masses m2 = mi/Ai, the 
pieces of both dispersions VLi{k) and VL2{k) which are proportional to k"^ vanish, as 
should be evident from (58) - (61). Thus, to order A;^, Q,2{k) vanishes identically. 



was stated above, here we assume that the coefficient of |fc| in the first equation in (58) is 
positive. In case it is negative, the term linear in would appear in Q.2{k) instead of f2i(fc). 
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while Qi{k) attains the value of the dispersion corresponding to the single-family 
Calogero model at its fermionic point A = 1. In other words, at this particular 
point in parameter space, our two-family Calogero model appears to be similar (but 
not entirely equivalent) to a system of free fermions. The vanishing of (A;) at 
/c 7^ is an indication of an instability at that point in parameter space. Evidently, 
one has to go beyond the quadratic approximation to resolve this problem. 

To summarize our discussion thus far, we have found that the low-momentum 
behavior of the two-family Calogero model is effectively that of an ordinary single- 
species Calogero model and a decoupled system of free massive bosons. This conclu- 
sion will manifest itself again when we compute density-density correlation functions 
in the next section. 

Let us discuss now how to determine the effective one-family parameters m, A and 
p. Eqs. (59) and (60) comprise two relations among the three effective parameters 
p, m and A . The third relation, needed to determine them unambiguously, arises from 
identifying the ground-state energy (31) with that of the single-species collective field 

namely, 

(AiPi,o + Ai2P2,o)' = -^P^ ■ (62) 



6Aimi ' ' 6m 

Note that although, strictly speaking, the first equality in the expression (31) for 
the ground state energy is based on the leading large- A'" behavior of the effective 
potential (27), the result quoted there is, in fact, the exac^ value of the ground state 
energy, valid also for finite A^i,2- Indeed, the integrand appearing in that equation 
(namely, the LHS of (62)) is invariant under the duality transformations of Section 



"That is,f7(fc) = i^fc2 + ^|fc|. 



2m 

^^Clearly, the free quasi-particles described by i}2{k) cannot contribute to the ground state 
energy. 
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3, which hold also beyond leading order. Thus, we can treat (62) on the same footing 
as (59) and (60), as far as large- A?" accuracy is concerned. 

Recall that the couplings and masses are related according to (2). Thus, elim- 
inating, for example, A2 and A 12 in terms of Ai and the ratio of masses, and then 
substituting these into (59), (60) and (62), it is possible to cast these three equa- 
tions into forms from which elimination of A and p in terms of m leads to the simple 
relations 

A = ^ , mp = mipi,o + m2P2,o , (63) 

m"^ mi 

while m itself satisfies the quadratic equation 

777 

+ ((1 - Ai)pi,o + (1 - A2)p2,o) — = . (64) 

(Note that Eqs.(63) and (64) are invariant under the duality transformations of 
Section 3, having been derived from duality-invariant equations.) 

After solving (64) for m we can obtain A from (63). It is natural to conjecture 
that the quasi-particles of the effective one-family Calogero model have fractional 
statistics with statistical parameter A. It will be interesting to prove this conjecture. 

As an example for the discussion above, consider the case in which all couplings 
and masses are positive. Under these conditions, (64) has two real solutions of 
opposite signs (which either remain invariant or are flipped by the duality trans- 
formations). Let us choose the positive root as desired mass m. Substituting the 
latter into (63) we obtain A and p, which are also positive. (We do not present the 
corresponding expressions here.) 

Finally, note that A = Ai, m = mi, subjected to the constraints (22) defining 
the SOHI, namely, the situation corresponding to (23), is also a solution of (63) and 
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1712 

Pl,0 + P2,0 

mi 



mi/ 



(64). 



6 Ground-state wave-functional and correlation 
functions 

Having completed diagonalization of Hcoii in (57), we next compute the correspond- 
ing collective- field ground-state wave functional. The ground state |0) is defined 

by bi{k)\0) = 62(^)10) = 0, for every k. Due to the orthogonal transformation 
(49), the latter conditions are equivalent to ai{k)\0) — a2{k)\0) — 0, or in real 
space, to 

Ai{x)\0) = A2ix)\0) = . (65) 

We shall now concentrate on the latter representation. Following the definitions (34) 
and (35) we write (65) as 

V dr]i{x) 2pi,o J x-y J x-yj 

V dr]2{x) 2p2,o J x-y J x-y J 

(66) 

where ^o[vi:V2\ — ({?7i, ?72}|0) is the wave functional we sought for. The solution of 
these linear and homogeneous coupled functional differential equations is evidently 
a Gaussian in the fluctuating flelds rji and r]2, which we readily flnd as 



x) - + Y In |a; - y\j r]i{y) 



+V2{x) \ - y) + Y In k - y| j miy) + Xi2r]i{x) in \x - y\r]2{y) \ , (67) 
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where the normaUzation constant Zq is fixed by the requirement that / T>if]iDif]2\^ ^[r]!., '/72]P 
1 . We may write (67) more compactly as 

2 



*ok,^2] = -^exp(-i / / dxdy ^ r]a{x)Kah{x - y)r]i{y)\ 



(68) 



where 



Ku{x) = ^5(a;) -2Ailn|a;| , 

Pi,o 

K22{x) = l:i^6{x) - 2\2\n\x\ , 

P2,0 

Ku{x) = ^21(2;) = -2Ai2ln|x|. (69) 

By comparing (36) and (69) we obtain the simple relation 

[Aa{x),Ai{y)] = d,dyKab{x - y) . (70) 

Now that we have the collective-field vacuum wave-functional, we can calculate the 
density-density correlation functions 

{Q\r]a{x)'qb{y)\Q) = ^/ T^ViT^m Va{x)rjb{y) exp(^-^ j j dxdy ^ 7]c{x)Kcd{x-y)r]d{y) 

(71) 

by standard quantum field theoretic methods. We thus obtain 

{0\Va{x)v,{ym = K^,\x - y) , (72) 
where K^,j^{x — y) denotes the inverse kernel 

/ dzJ2Kac{x - z)K^f^\z -y) = SabS{x - z) . (73) 

^ c=l 

Throughout the discussion from (67) to (72) we have assumed, of course, that 
-f^a6(a;—y) is positive definite. Since the kernels Kab{x—y) are translation-invariant, 
it is more appropriate to calculate the inverse kernels in momentum space. Thus, 
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by Fourier transforming both sides of (70), we find, using (38), (40) and (43) - (47) 
that 

/oo 
dke^'^Kabik), (74) 

with 

Ka,{k) = J^^'^- (75) 
Note that Kab{k) diverges hke l/\k\ as A; — > 0. Eq. (75) can be also verified by 
Fourier transforming all equations in (69), by employing the Fourier integral 



rria nib ^ab{k) 



In \x\ 



ikx 



k\ 



(76) 



IJ — oo 

(and the standard Fourier representation of 5{x)). The desired inverse kernel is then 
found to be 



K-,\x-y)^ 



1 



dke'^^^-y^K-,\k), 



(77) 



where K ^{k) is the matrix inverse to (75). ^From (75) and (44) - (47) we obtain 
the latter as 



27rpi,()P2,0 



V 



-A 



12 



1-Ai 



2np 



1,0 I I ^ / 



where 



K = (1 - Ai)A2P2,o + (1 - A2)Aipi,o . 
In obtaining (78) we used 

;i-Ai)(i-A2) 



detK(A;) = 



1 

27rpi,oP2,o|^| 



K + 



271 



k\ 



(78) 



(79) 



(80) 



where (2) was invoked. By substituting (78) in (77) we obtain the integral repre- 
sentation 



{0\rjaix)vbiOm-Kbi^)- 



Pi.o 
1-Ai 



27r 







P2,0 
I-A2 



5{x) + 
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- 2 



-^12Pl,0p2,0 \ 

(i-Ai)(i-A2) y 
/ dk 

Al2Pl,0P2,0 -^2p2_o 

V (1-Ai)(l-A2) 



cos{kx) 



k + 



2-KK 



(81) 



(1-A2)2 



(1-Ai)(l-A2) 



for the desired correlators (72). The latter integral is well defined for 

27r« 2j:^ + :^Uo, 



(82) 



" (1-Ai)(l-A2) - VI -Ai ■ l-\2, 
in which case it can be expressed as a combination of the Sine-integral si(2;) = 

oo oo 

— / dt{smt/t) and the Cosine-integral ci(2;) = — / dt{cost/t) [58], leading to the 

2 Z 

final expression (for x ^ 0) 

( 

(0|77„(x)776(0)|0) = 2 



Al2Pl.(]P2,0 \ 

(1-Ai)(l-A2) 



Al2Pl,0P2,0 

V (1-Ai)(l-A2) 



^2P2,0 
(1-A2)2 



F{q\x\ 



(83) 



where 



(84) 



F{q\x\) — sin {q\x\) si {q\x\) -\- cos {q\x\) ci {q\x\ 

Eq. (83) generalizes a similar result known for the one- family Calogero model [56] . 
The function F{q\x\) has oscillatory behavior for small to moderate values of its 
argument (as one would typically expect of correlation functions), and falls like 
l/{qxY at large values. In fact, F[q\x\) has a rather simple asymptotic expansion 
[58] in inverse powers of (gx)^: 

(-l)"(2n- 1)! 



n=l 



{qx) 



2n 



q\x\ » 1 



(85) 



In the case q < Q the integral representation (81) is ill-defined, due to the pole 
on the integration contour, however, the asymptotic series (85) still makes sense. 



^■^It is straightforward to obtain this asymptotic expansion by repeatedly integrating-by-parts 
the integral in (81). 
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as the large- 1 a; I behavior of the correlation function comes from momentum values 
I A; I << In either case, we can trust the asymptotic series 

{0\rja{x)rj,mO) = K-,\x) 



(1-Ai)2 (i-Ai)(l-A2) 



. 2Al2pl,0p2,0 ^'^2^2.1) . 

\ (1-Ai)(l-A2) (1-A2)^ / 

2t:k\x\ 



E(-ir(2.-i) f^^^'Jl^-^-' ]".(86) 



n=l 



valid for (i__x^)(i-a2) -'- ' capture the leading long-distance behavior of the 
correlators. 

Close inspection of (83) reveals that the coefficient matrix there is of rank-1. In- 
deed, due to (2) it can be written as ww^ with w„ = . Thus, det (0|77„(x)776(0)|0) 
for all a; 7^ .^"^ More importantly, this means that only one combination of the 
density fluctuation fields 77a (a;) can have long-range correlations. In order to reveal 
the reason for this behavior, it is instructive to rewrite (68) in momentum space. 
Starting from the identity 

^ = j j ^^^y J2Va{x)Kabix - y)rihiy) = i'^'^f j dk Ylfjai-k)kab{k) fjb{k) 

a,b a,b 

we obtain, using (75), (50) and Uab{—k) = Uab{k) (see (54)) that 



S^{27^fj dkj:U-k)-^Uk), (87) 



where 



Uk) = j:uab{k)J^fjb{k). (88) 

T V ^Pbfi 

Thus, (68) factors into two uncorrelated Gaussian pieces, describing the vacuum 
fluctuations of the two fields Ca{x) = J dke^^^(a{k)- ^From the low-momentum 



Related to this is the fact that the matrix K ^(0) exists, but due to (2) it is not invertible. 
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behavior (58) of Qi{k) and (/c), and from (59) - (61), we conclude immediately that 
at long-distances, the vacuum fluctuations of the combination Ci(^) coincide with 
those of a single-species Calogero model with parameters A, m and p, whereas C2ix) 
appears like a white-noise random field with ultra-local correlations proportional to 
S{x — y). All this, of course, is in full accordance with our earlier conclusion that the 
two-family Calogero model is equivalent, in the large- limit and at long-distances, 
to two decoupled systems, one - an effective one-species Calogero model, and the 
other, a system of free massive particles. 

As a simple application of this observation, consider the special two-family sys- 
tems living on the SOHI (22) in parameter space. Let us use our result (83) to 
evaluate the correlation function of the effective one-family density (23). We shall 
denote the fluctuating part of p{x) in (23) by ri{x) = r)i{x) — jr]2{x). Substituting 
the parameters corresponding to (22) and (23) in the appropriate places, we obtain 

(0|r^(a;)r^(0)|0) = (0| (r^i(x) - jrj2{x)^ (r^i(O) - ^7^2(0)) |0) 



{0\Vi{x)vimO) - j{0\vi{x)v2m0) + ^(0|7?2(x)772(0)|0) 



2A / 1 \ ^ 
YZ^[Pi,o-jP2fi) F{q\x\), (89) 



with 



2^^ / 1 \ 

Q = ^— Y yPifi - jp2,o J ■ (90) 

^From (85) we obtain the leading long-distance asymptotic behavior of this correla- 
tor simply as 

(0|r/(a;)77(0)|0) = -^^, (91) 
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which does not depend on the particle density at all, in accordance with the one- 
family correlation function of [57, 59]. 



6.1 The low- lying collective excitations 

We close this section by analyzing of the low-lying collective excitations above the 
vacuum. These excitations are the pseudo-particles created by bi 2{k) acting on the 

vacuum. The corresponding states are 

1 °° I 

biim = ir^ I dxe^'^j:ua,{k)J^Al{xm, (92) 
2^/71 J ^ \ nib 

where we used (49), (38), (34), (35) and the fact that the Uab are real. Next, we 

observe from (34), (35) and (67) that 

4(,)|0) = -2, ( + A J**M + A. J*5£M) 10) (93) 
\ 2 pao J x-y J x-y J 

where a is the index obtained from a by permutation, i.e., 1 = 2,2 = 1. By 
substituting (93) in (92), and using (39), we then obtain, after some algebra, 

bl{k)\0) = ^ / dxe^'^Y.^ab f./^u;,5 77,(x) + MubiVlix)] |0) , (94) 
J-nk J V V \ Pm \ Phn 



PbO V PbO 

where we also used (44) - (46) on the way. Performing summation over the index 
b and using (50), which is uquivalent to the relations Uaa^^aa + ^aai^aa — ^a^aa and 
UaaOOaa + UaaOJaa = ^aUaa , we may simphfy (94) further into 



oo 



biikm^'-^ f dxe^'^{J^u^^r)^{x) + J^u^,ri,{x)]\0). (95) 

— oo ' ' 

Using the definitions of Uab given above (55) we thus obtain 



bUkm oc ,/^^^i^cot0 / dxe'''''r]i(x)\0)+ [ dx e'''''r]2(x)\0) 
\l m2Pi,o J J 



bl{k)\0) (X /^^^i^ tan</) / dx e'''''r]i(x)\0) + [ dx e'''''r]2{x)\0) . (96) 
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Let us now compute cot 0. Substituting (55) in the identity cot = cot 20 + 



\Jl-\- cot^ 20 and then using (44) - (46) we obtain 

cot = , (97) 

which upon expansion to leading order in yields 

' \ mi mi I 

Finally, returning to the original configuration space coordinates 

^e'''^^ , / dxe'''''rj2{x) ^ 

=1 •' a=l 

and substituting it and (98) in (96) we find 



/ dxe'''''rji{x) ^e'''^^ , / dxe'''''rj2{x) ^ e'''^" , (99) 



„, / 1 1-^1 _ 1-^2 \ Ni N2 



6l(A;)|0) oc ^ 1 + Ee^'"|0) + E e''^"|0) (100) 

\ mi m2 / 

and 



hi{km oc 1 - ^ Co^Ar;.o i^i E^^'^io) + e ^^'^io) , (loi) 

where we used A1/A12 = mi/m2 from (2). Our low-energy excitations are in fact 
long- wavelength density oscillations (phonons) in which the velocities of all particles 
point in the same direction, either positive or negative, depending on the sign of k. 

Comparing the excitation (100) with that obtained in [36], we see that our rela- 
tive weight approaches the relative weight m\/m2 from [36], but only in the small 
k hmit. The second excitation (101), which describes the free-particle motion, does 
not have any counterpart in Sen's approach 
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7 The multi-species generalization 



Let us now apply the methods we used in our analysis of the two-family Calogero 
model to investigate its F-family generahzation. Thus, consider a collection of F 
species of particles. The a— th family contains Na particles of mass rUa (a = 
1, 2, . . . F) , which interact among themselves with coupling constant = Aa^. The 
inter-family mutual interaction strengths are \ab — Ma 7^ ^) ■ All couplings 
satisfy the constraint 

^"^ (102) 



mama 

where the same constant c holds for all pairs of indices. This condition generalizes 
(2) and guarantees that there be no three-body interactions [43, 44]. It follows from 
(102) that 

Kh = \lb, (103) 

similarly to (3). We take the limit in which all particle numbers Na tend to infinity 
at the same rate. The collective-field Hamiltonian of the F— family Calogero model 
is given by 

Hcoll = TT— I dx \d:,'Ka{x)pa{x)d:,'Ka{x)+pa{x) i ^"' ^ + ^ Xabi ^^^^^^ 

a=l \ \ Pa b=l ^ y , 

(104) 

up to singular boundary-contributions (compare to (16)). The collective density 
fields and their conjugate momenta obey canonical equal-time commutation rela- 
tions analogous to (12). The a— th density field is subjected, of course, to the 
normalization condition / dxpa{x) = Na- 

The ground-state energy and density distributions can be found, as in the 2- 
family case, by varying the effective potential corresponding to (104) with respect 
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to all Pa ■ We thus obtain the uniform ground-state densities 

N 

Pafi^^. (105) 

where L is the length of the large confining box, and the ground-state energy 



-Tr2^2 / F 

hlQ — 



6L2 

In the last equation, which is a direct generalization of (31), we used (102). 

Similarly to our conventions in the F — 2 case, we shall take all Aa > 0, 
and anticipating the possibility of negative masses and densities, we shall also set 
sign nia — sign o for all families. 

We now expand around the uniform configurations (105) 

pa{x) = pa,0 + r]a{x), (107) 

and introduce the F operators Ac{x) analogous to (34) and (35). The Fourier 
transforms adk) of these operators, and their hermtian conjugates, defined as in 
(40) and (43), obey the commutation relations [aa{k) , a\{k')] — uJab{k)S{k — k') , 
where the F x F real symmetric matrix 



U^a,(k) = ^—^kX, + Xa,7r.h^\k\ (108) 

2ma V rriamb 

is an immediate generalization of (47). 

Let us now diagonalize ujab and obtain its eigenvalues fla{k). We shall contend 
ourselves with the low-momentum behavior of these dispersions, and compute them 
to order fc^. Details of this computation are given in Appendix B. We find that one 
of these eigenvalues is given approximately by 

"■W- ;V y + '^E^W- (109) 

^i=a ma «=1 " 
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Since in our conventions Xa > and signm^ = signpa^Oi the coefficient of in the 
last term is manifestly positive. The remaining F — 1 eigenvalues fl2{k) ~ ^3{k) ~ 
. . . ~ Qpik) ~ k'^. Thus, as in the F — 2 case, the first dispersion 0,i(k) looks 
effectively like that of a one- family Calogero model. The remaining F—1 dispersions 
correspond evidently to free massive bosonic quasi-particles, decoupled from each 
other and from the Calogero-like mode. 

Let us now compare our Calogero-like dispersion law (109) with that given by 
Eq.(58) of [36], which was derived to leading order in j. The terms linear in 
momentum k agree, as can be easily seen by virtue of the relation 

, (110) 

where vp denotes the "Fermi" velocity in Sen's approach. We observe that the 
k^ correction in our dispersion (109) arises from both the linear and the quadratic 
differential operators in the Hamiltonian (104). Taken together, they give a con- 
tribution of the typical weight ^^y^. These terms have been completely ignored in 
Sen's corresponding dispersion relation mentioned above, because they did not seem 
to be analytically computable in that approach. In this respect, our method seems 
to do better than the generahzed Thomas- Fermi method of [36] . 

7.1 Duality transformations and symmetries of the F- family 
model 

The duality symmetries of the two-family model, discussed in Section 3, have a 
natural generalization in the F-family model. The main idea in searching for these 
duality symmetries in the generic case is to find all classes of physically equivalent 
F— family models. Here, the eqivalence of two systems in a given class means that 
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despite the fact that these systems are defined in terms of different sets of collective 
fields, momenta and parameters, their collective Hamiltonians (104) coincide. 

As our first step in classifying all these duality symmetries, it is easy to check 
that the Hamiltonian (104) is invariant under the following set of elementary duality 
transformations, associated with a prescribed family index a: 

Xa = Xab = for 6 7^ o; Xcd = Xcd: for c,dy^a; 

rua ^ „ , , 

rria = — r— ; = m^, for 6 / a; 

Pa = -KPa, Pb = Pb, for b^a; 

IT 

TTa = -T^; TTft = TTfe, for 6 7^ a. (Ill) 

Let us denote all these transformations collectively as Ta , a = 1, . . . F . As can be 
readily non-trivially only on the a— th family. Clearly, = I , where / 

is the identity transformation. Therefore, and / form a Z2 group of of symmetries 
of (104). It is possible to compose any pair of such elementary transformations, 
corresponding to different family indices. Moreover, this composition is obviously 
commutative. In this way we can generate the more complicated duality symmetries 
TaTb — TbTa of (104). Continuing this process, we can generate all transformations 
involving triplets of different family indices TaT^Tc all the way up to the maximally 
possible composition T1T2 ■ ■ - Tp. In this way, the elementary transformations T^, 
together with the identity / generate an Abelian group of order 2^ , which exhaust 
all possible duality symmetries of (104). Evidently, each element in this group is of 
order 2. We can readily identify this group is as ®^ Z:2 ■ 

Within this general picture, we can interpret the two-family duality transforma- 
tions (20), (21) as an elementary transformation on the second family, while the 
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transformations (18), (19) may be understood as the simultaneous action of the two 
elementary transformations on the first and the second family. 

All physical properties of the F-family system should remain invariant under 
the duality symmetries of (104). In particular, the F dispersions ila(^) should be 
duality-invariants. This can be proved in a manner similar to the F — 2 case, 
namely by showing that the effect of a duality transformation on the matrix uj in 
(108) is simply a conjugation by an orthogonal matrix. It is enough to prove 
this conjugation for an elementary transformation (111): It can be easily checked 
that the elementary duality transformation leaves u!aa,'^bb and cUcd {i>,c,d,^ a) 
invariant, while cUab — ^ba Aips sign. All this can be summarized by saying that 
under Ta, cu ^ Ta^^a, where the diagonal orthogonal matrix is obtained from the 
unit matrix by flipping the sign of its a-th diagonal entry. Thus, all the dispersions 
fla{k) are invariant under the duality transformations of (104). In particular, all 
duality transformations preserve positivity of the coefficient of in the second 
term in (109). 

7.1.1 Special F-family Calogero models which resemble single-family 
Calogero models 

In Section 3.1 we applied the duality group to identify a special subset of two- 
family models, with apparently different families, which nevertheless resemble, in 
almost all respects, a Calogero system composed of identical particles. Those are 
the two-family systems living on the so-called "surface of hidden identity" (SOHI) 
(22) in parameter space. Not surprisingly, this construction can be generalized to 
the F-family case, as we now show. 
^^See the discussions following (47) and (51). 
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In order to find special multi-family Calogero models, which resemble a single- 
species Calogero model with some effective particle density p, mass m and 
interaction strength A , let us spht the latter effective one-family system into F 
different sub-families whose particles share the same mass m, the same coupling 
constant A and the same mutual interaction strength Aa(, = A. This splitting means 
of course that we no longer consider totally symmetric many-body wave functions, 
and replace them by wave functions symmetric separately in the coordinates of 
the particles of each family. Let pa denote the density of a-th family particles. 
Therefore, p — J2aPa- We shall refer to the F-family system thus obtained as 
maximally degenerate (MD). 

Let us consider our MD F-family model as the image under one of the elementary 
duality transformations in (111) for some 1 < k < F. Since T| = /, its origin 
must have been the less-degenerate F-family Calogero system with parameters 

Ajfc = — , Xak — Afea = — 1 , Ao = Aab = A for a,b ^ k; 

in 

rrik — — — , rUa — m for a^k; (112) 
A 

partial densities 

-y, Pa for a^k, (113) 

and total density 

p{x) = E - \pki^) ■ (114) 

To reverse the argument, the k-th family in the less degenerate F-family system de- 
fined by (112) - (113) is singled out, yet the latter is dual to the MD F-family system, 
and therefore resembles a single family Calogero model. We can now apply another 
element Ti (/ ^ k) of (111) on the less degenerate F-family system, and obtain an 
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even lesser-degenerate F-family system, now with two families singled out. Yet, the 
latter system is stil dual to the MD system. We can clearly continue this process 
thus generating all images of the MD F-family system under all 2^ — 1 non-unit 
elemenets of the duality group. Each of these images comprises a special F-family 
system, whose families are not all identical, which nevertheless resembles a single 
family Calogero model. Just to give a nontrivial example, consider the image of the 
MD F-family system under the combined duality transformation Tfc_|_iTfc_|_2 ■ ■ - Tp. 
This image is the F-family Calogero system with masses m^, coupling constants 
Aa and mutual interaction strengths given by 

Ai = A2 = ... — \k — \] 

mi = m2 — ... = ruk — m; 

Xab = A, for a,b <k; Xk+i = Xk+2 = ■■■ = Af = 

m 

vfik+i = mk+2 = ■■■ = mF = -—; 

^ab — T") for a,b > k; Xab — for a < k and b > k. (US) 
A 

The effective one-family collective-field is obviously given by 

k I ^ 

a=l a=k+l 

For F — 2 this result clearly coincides with (23). 

8 Concluding remarks 

We conclude this paper with a few comments on some recent papers and also on 
another set of duality transformations which were discussed in the literature. 
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Sergeev and Veselov [39] constructed supersymmetric extensions of the Calogero- 
Sutherland model which actually correspond to our two-family Calogero model with 
A1A2 = 1, A12 = —1 and m\m2 < 0! They gave solutions in terms of deformed Jack 
polynomials. In a recent paper, Kohler and Guhr [41] introduced a supersymmetric 
generalization of the Calogero-Sutherland model. Their construction is based on 
Jacobians for the radial coordinates on certain superspaces. This approach allowed 
them to explicitly construct the solutions in terms of recursion formulae for a non- 
trivial (A1A2 = 1) one-parameter subspace in the (Ai,A2) plane. The underlying 
model involves two kinds of interacting particles, one with the positive and the 
other one with the negative mass. Needless to say, this again corresponds to our 
two- family Calogero model with A12 — —1. It is interesting to observe that the 
authors of Refs.[39]-[42] were probably unaware of the constraints (2). Namely, in 
their approaches, these constraints remain hidden, but still present, as can be easily 
checked by direct substitutions. Consequently, the two types of models discussed in 
[39]- [42], share the very same parametric structure, which enables one to transform 
them to the one-family Calogero model. This connection then guarantees their exact 
integrability. Although our collective-field approach is applicable only to the multi- 
species Calogero system with an infinitely large number of particles within each 
family, we believe that our findings shed some light on the problem of their exact 
integrability in general. 

Let us finally remark on an additional approximate duality invariance of the 
multi-species model (104), known in the literature [60, 61]. For simplicity, we illus- 
trate it on the two-family Hamiltonian (16), in which case this approximate duality 
symmetry involves special two- family systems, in which A2 = ij^ but A12 = +1, as 
opposed to (22). This duality symmetry is destroyed by subleading terms in (16) as 
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we now explain. 

Estimating the dependence of the terms in the effective potential (25), we 
see that the dxlnp\ and dxlnp2 terms are down by and by respectively, 
when compared with the Hilbert-transform terms. The effective potential to this 
order is then 

zmi J \ J X — y J X — y J 



^ 2 

21712 J \ J X — y J X — y 

One can verify that the above potential is invariant under the following set of trans- 



formations of the parameters: 

A2 A2 ^^2 

and of the densities 

Pi = P2; P2 = ±A2P2- (118) 

Arguing along the same lines as before, we can show that the two-family Calogero 
model is equivalent to the one-family Calogero model if the following conditions are 
satisfied: 

1 ni 
M = A; A2 = -; mi = m, m2 = ±—] \i2 = ±1, 

where A denotes the coupling constant and m the mass of the effective one-family 
model. The effective one-family collective density is given by 

p{x)=pi{x)±jp2{x). (119) 

We stress that the above picture is valid in the leading approximation only. Conse- 
quently, this description is inevitably destroyed by the next-to-leading-order terms 
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stemming, for instance, from the quantum collective-field fluctuations. In contrast, 
the various dualities Ti, T12 (and, of course, the indentity /), discussed in Section 
3, are exact symmetries of the collective Hamiltonian (16). Consequently, the dual 
equivalence of any two-family system on the SOHI (22) to a two-family system of 
identical families given by (23) or (24) is also an exact property of (16). 
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Appendix A: An alternative derivation of (27) 



For completeness, we present here an alternative derivation of the expression (27) for 
the collective effective potential. To this end, it is useful to introduce the wieghted 
density 

p{x) ^ ^fx^pi{x) + ^fx^P2{x) , (A.l) 

and consider the resolvent 

Giz)^[^ (A.2) 
J z-y 

associated with it, in which z is a complex variable. As in the single family case, the 
particles are expected to condense in the ground state in a large but finite segment 
of some length L - the support of p{x). The resolvent G{z) is therefore analytic in 
the complex plane, save for a cut along that segment, which lies on the real axis. 
Consequently, as usual, we have 

G{x T iO) = F{x) ± mp{x) (A.3) 

where 



J x-y 



y 

The combinations of principal part integrals in the first two lines in (25) are evidently 
proportional to F{x). This fact, together with (3) can be used to rewrite (25) as 

8mi J pi 8m2 J p2 



+ / dx F{x) — - — pi H — p2 



Till J V 2 



+ 1^1 (iVi - J dxpi{x)j + p2 (^N2 - J dxp2{x)j . (A.5) 
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The term p{x) F'^{x) in the third hne in (A. 5) appears to be a trilocal cubic functional 
of the densities. In fact, it can be brought to local form by a standard trick [51], 
based on the analytic structure of G{z) and its asymptotic behavior 

G(.) ~ ^■^ + ^^^ (A.6) 

z 

as 2; — > oo. Note that the fact that pi{x) and p2ix) have compact support is 
necessary for (A.6) to hold. It follows from (A.6) that dzG^{z) = 0, where C^o 
is a circle around the point at infinity. Then, shrink to wrap around the cut to 
obtain 2i J dxQG^{x — iO) — 0, namely, 

2 

I dxp{x)F''{x) = y / dxp^x) . (A.7) 
Substituting (A.7) in (A. 5) we obtain the effective potential in its final form as 



-P2 



<,mi J pi »m2 J P2 

^ f rn/__X o - 1 , A2 - 1 

mi 



+ ^ fdxF{x)dJ^^p,+ 
+ 1^1 (a^i - / dxpi{x)^ + p2 (a^2 - J dxp2{x)^ , {A.t 



which can be shown easily to coincide with (27). 
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Appendix B: Diagonalization of uJah{^) in (108) 



The matrix ujah{k) in (108) is of the general form 

u; = DA;2 + vv'^|A;|, (B.l) 
where D is a diagonal matrix with diagonal entries da — and v is an F- 



dimensional vector with components Va — cniayj^^^- — sign (ma) y^^^^a! ^ where 
we used (102) twice. After a standard computation we obtain the characteristic 
polynomial of (B.l) as 

Pf{z) = det(^ - D/t' + vv^|A;|) = det(^ - DA;') ( 1 - — ^Ir^v ) . (B.2) 

In the generic case, where all the diagonal elements da of D are different from each 
other (and we shall henceforth assume this case), none of the eigenvalues of uj is an 
eigenvalue of DA;'. Therefore, all roots of Pf{z) are obtained by finding the zeros 
of the second factor in (B.2), namely, from 



T 

V 



;v = E-^ = i. (B.3) 



z — T>k'^ ~^ z — dak"^ 

It is clear from (B.l) that all eigenvalues of uj vanish at least as fast as |A;| when 
the latter tends to zero. Alternatively, this can be seen by taking the hmit |A;| — > 
in (B.3). Since all v1 > 0, the roots must vanish at least hke |A;|. Thus, let us set 
^ = C|A;| in (B.3), and solve 



^®This should be well-known. The reason for this is simple: the second factor in (B.2) has 
simple poles located at the eigenvalues of Dfc^, which cancel against zeros at the same points in 
the first factor in (B.2). If D is not degenerate, all those zeros are simple, and we end up with a 
nonvanishing result. 
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for (. In fact, we need to find ( only to first order in \k\, as we are interested in the 
eigenvalues of (B.l) only to order k"^. 

Let us first look for roots of (B.4) which do not vanish as —>■ 0. We find only 
one such root, which at = is given simply by = Ea'^l > 0- To find the 
leading correction to this root we expand (B.4) in inverse powers of ( to order |A;|: 

l-7E^a-||E^X + ^(^')=0. (B.5) 
^ a=l ^ a=l 

Only one root of this approximate quadratic equation converges to , and to order 
I A; I it is given by 

^(1) = j2vl + \k\ . (B.6) 

a=l ^a=l '^a 

Upon substituting the appropriate values of Va and da, (B.6) coincides with the 
expression (109) quoted in the text. In our conventions Aq > and sign ma = 
sign pafl . Therefore all Va are real, and consequently the coefficient J2a=i '^a 1^1 
in (109) is positive. 

All other F — \ roots of (B.4) vanish as \k\ when the latter tends to zero. 
To see this, restore the original roots z = C,\k\ and consider det(— a;) = -Pf(O) ~ 
+ \ k\v'^'D ^v^ . As — > this behaves like (k'^)^ ^ m ■> proving our 

claim. 
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